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Abstract.  We  prove  new  velocity  averaging  results  for  second-order  multidimensional  equa¬ 
tions  of  the  general  form,  C(S7X,  v)f(x,  v )  =  g(x,  v )  where  C(X7x,v)  :=  a(v)  ■  Vx  —  Vj  ■  b(v)Vx. 
These  results  quantify  the  Sobolev  regularity  of  the  averages,  f  f(x,  v)<f>(v)dv,  in  terms  of 
the  non-degeneracy  of  the  set  {u  :  |£(i£,  u)|  <  <5}  and  the  mere  integrability  of  the  data, 
(/iff)  £  (^x  vi  Lx  v)-  Velocity  averaging  is  then  used  to  study  the  regularizing  effect  in  quasi- 
linear  second-order  equations,  C(ff7x,p)p  =  S{p)  using  their  underlying  kinetic  formulations, 
C(Vx,v)Xp  =  ffs-  In  particular,  we  improve  previous  regularity  statements  for  nonlinear  con¬ 
servation  laws,  and  we  derive  completely  new  regularity  results  for  convection-diffusion  and 
elliptic  equations  driven  by  degenerate,  non-isotropic  diffusion. 
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1.  Introduction 

We  study  the  regularity  of  solutions  to  multidimensional  quasilinear  scalar  equations  of  the 
form 

d  f)  d  r)2 

(L1)  = 

j=0  J  j,k=  1  J 

where  p  :  M1+d  — »  M  is  the  unknown  held  and  Aj,  Bjk,  S  are  given  functions  from  M  to  M. 

This  class  of  equations  governs  time-dependent  solutions,  p(t  =  x0,x),  of  nonlinear  conser¬ 
vation  laws  where  A0(p)  =  p  and  Bjk  =  0,  time-dependent  solutions  of  degenerate  diffusion 
and  convection-diffusion  equations  where  {Bjk}  >  0,  and  spatial  solutions,  p(x),  of  degenerate 
elliptic  equations  where  Aj  =  0.  The  notion  of  solution  should  be  interpreted  here  in  an  appro¬ 
priate  weak  sense,  since  we  focus  our  attention  on  the  degenerate  diffusion  case,  which  is  too 
weak  to  enforce  the  smoothness  required  for  a  notion  of  a  strong  solution.  Instead,  a  common 
feature  of  such  problems  is  the  (limited)  regularity  of  their  solutions,  which  is  dictated  by  the 
nonlinearity  of  the  governing  equations.  A  prototype  example  is  provided  by  discontinuous 
solutions  of  nonlinear  conservation  laws.  In  [LPT94a],  Lions,  Perthame  &  Tadmor  have  shown 
that  entropy  solutions  of  such  laws  admit  a  regularizing  effect  of  a  fractional  order,  dictated 
by  the  order  of  non-degeneracy  of  the  equations.  In  this  paper  we  extend  this  result  for  the 
general  class  of  second-order  equations  (1.1).  In  particular,  we  improve  the  [LPT94a]-regularity 
statement  for  nonlinear  conservation  laws,  and  we  derive  completely  new  regularity  results  for 
convection-diffusion  and  elliptic  equations  driven  by  degenerate,  non-isotropic  diffusion. 

The  derivation  of  these  regularity  results  employs  a  kinetic  formulation  of  (1.1).  To  describe 
this  formulation  let  us  proceed  formally,  seeking  an  equation  which  governs  the  indicator  func¬ 
tion,  Xp{x)(v)  '■=  sgn(u)(|p| —  M)+  associated  with  p,  and  which  depends  on  an  auxiliary  velocity 
variable  v  E  M,  borrowing  the  terminology  from  the  classical  kinetic  framework.  To  this  end, 
we  consider  the  distribution  g  =  g(x,  v),  defined  via  its  velocity  derivative  dvg  using  the  formula 

(1.2)  dvg(x,v)  :=  (a(u)  •  V*  -  Vj  •  b(v)Vx  +  S(v)dv^xP(x)(v),  a j  :=  A'-,  bjk  :=  B'jk  >  0. 

Observe  that  the  nonlinear  quantities  <f>(p)  can  be  expressed  as  the  u- moments  of  yp,  <f>(p)  = 
f  $'(v)xP(v)dv,  $(0)  =  0.  Therefore,  by  velocity  averaging  of  (1.2)  we  recover  (1.1).  Moreover, 
for  a  proper  notion  of  weak  solution  p,  one  augments  (1.1)  with  additional  conditions  on  the 
behavior  of  <f>(p)  for  a  large  enough  family  of  entropies  <l>’s.  These  additional  entropy  conditions 
imply  that  g  is  in  fact  a  positive  distribution,  g  =  m  E  A4+,  measuring  the  entropy  dissipation 
of  the  nonlinear  equation.  We  arrive  at  the  kinetic  formulation  of  (1.1) 


(1.3)  £(Vx,v)f(x,v)  =  dvm  -  S{v)dvxP (x)(v),  f  =  XP ,  m  E  M+ 

where  C  is  identified  with  the  linear  symbol  C(if,v)  :=  a(u)  •  if  +  (b(u)£,£).  We  recall 
that  p(x)  itself  can  be  recovered  by  velocity  averaging  of  f(x,v)  =  Xp(u)>  yia  the  identity 
p(x)  =  f  f(x,v)  dv.  In  Section  2  we  discuss  the  regularity  gained  by  such  velocity  averaging. 
There  is  a  relatively  short  yet  intense  history  of  such  regularity  results,  commonly  known  as 
‘velocity  averaging  lemmata’.  We  mention  the  early  works  of  [GLPS88],  [DLM91]  and  their 
applications,  in  the  context  of  nonlinear  conservation  laws,  in  [LPT94a,  LPT94b,  LPS96,  JP02]; 
a  detailed  list  of  references  can  be  found  in  [Pe02]  and  is  revisited  in  Section  2  below.  Almost 
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all  previous  averaging  results  dealing  with  (1.3)  were  restricted  to  first-order  transport  equa¬ 
tion,  deg  £(*£,•)  =  1.  In  Section  2.2  we  present  an  extension  to  general  symbols  £’s  which 
satisfy  the  so  called  truncation  property,  luckily  —  as  shown  in  Section  2.4  below,  all  C s  with 
deg£(i£,  •)  <  2  satisfy  this  property.  If  the  problem  is  non-degenerate  in  the  sense  that* 


then  by  velocity  averaging  of  the  kinetic  solution  /  in  (1.3),  p(x)  =  fv  f(x,v)dv  has  a  Ws,r- 
regularity  of  order  s  <  (3a/ (3a  +  2)  with  an  appropriate  r  =  rs  >  1;  consult  (2.25)  below.  In 
the  particular  case  of  non-degenerate  homogeneous  symbols  of  order  k , 

3a  G  (0,1)  s.t.  sup  |{v  :  |£(i£,v)|  <  5}  I  <  Sa, 
kl=i 

we  may  take  (3  =  k  and  velocity  averaging  implies  a  -regularity  exponent  s  <  ka/ (3a  +  2). 
The  main  results  are  summarized  in  averaging  lemmata  2.2  and  2.3. 

In  Section  3  we  turn  to  the  first  application  of  these  averaging  results  in  the  context  of 
nonlinear  conservation  laws,  pt  +  V  •  A(p)  =  0,  subject  to  L°°-initial  data  po-  If  the  equation  is 
non-degenerate  of  order  a  in  the  sense  that 

|  {v  :  | t  +  a (v)  •  £|  <  h}  |  <  5a  and  sup  |a'(u)  •  £|  <  5l~a ,  V  r2  +  |£|2  =  1, 

(t>:  |r+a(u)-£|<5} 

then  for  t  >  0,  the  entropy  solution  p{t,x )  gains  Sobolev  regularity  pit,  •)  G  ITj(x!  (M/)  of  order 
s  <  sa  =  a/(2a  +  l).  This  improves  the  Sobolev- regularity  exponent  of  order  a/(a  +  2)  derived 
at  [LPT94a]  (while  facing  the  same  barrier  of  .s'i  =  1/3  discussed  in  [DLW05]). 

Section  4  is  devoted  to  convection- diffusion  equations.  We  begin,  in  Section  4.1,  with  second- 
order  degenerate  diffusion  pt  —  Yh^xjXkBjk{p)  —  0-  The  emphasis  here  is  on  non-isotropic 
diffusion,  beyond  the  prototype  case  of  the  porous  medium  equation  (which  corresponds  to 
the  case  when  Bj *.  is  a  scalar  multiple  of  the  identity,  Bjk  =  B5jk).  The  regularizing  effect  is 
determined  by  the  smallest  non-zero  eigenvalue  X(v)  =  A(b(u))  ^0  of  b(u)  :=  B'(v ),  so  that 

(1.4)  |  {v  :  0  <  X(v)  <  h}  |  <  Sa  and  sup  |(b,(u)^,^)|  <  h1-a,  V  |^|  =  1. 

{v:  A(t;)<(5} 

Staring  with  initial  conditions  po  G  L°°,  then  the  corresponding  kinetic  solution  p(t  >  0,  •)  gains 
1T/)X. -regularity  of  order  s  <  2a/(2a  +  1).  In  Section  4.2  we  take  into  account  the  additional 
effect  of  nonlinear  convection.  The  resulting  convection-diffusion  equations,  coupling  degenerate 
and  possibly  non-isotropic  diffusion  with  non-convex  convection  governing  capillarity  effects, 
are  found  in  a  variety  of  applications.  Consider  the  prototype  one- dimensional  case 

P,  +  A(p)x  -  B(p)xx  =  0.  A'(p)  ~  p’,  B'(p)  ~  |p|“ 

The  regularizing  effect  is  dictated  by  the  strength  of  the  degenerate  diffusion  vs.  the  convective 
degeneracy.  If  n  <  t  we  find  -regularity  of  order  s  <  2/(n  +  2),  which  is  the  same  Sobolev- 
regularity  exponent  we  find  with  the  ‘purely  diffusive’  porous- medium  equation,  i.e. ,  when 
A  =  0.  On  the  other  hand,  if  the  diffusion  is  too  weak  so  that  n  >  21,  we  then  conclude  with  a 
Sobolev-regularity  exponent  of  order  s  <  l/(£+2),  which  is  dominated  by  the  convective  part  of 
the  equation.  In  Section  4.2  we  present  similar  results  for  multidimensional  convection-diffusion 

*We  use  X  <  Y  to  denote  the  estimate  X  <  CY  where  C  is  a  constant  which  can  depend  on  exponents  such 
as  a,  /3,p  and  on  symbols  such  as  £,  Aj,  Bjk  but  is  independent  of  fields  such  as  p,  co-ordinates  such  as  x,  t,  v, 
and  scale  parameters  such  as  <5.  We  use  X  ~  Y  to  denote  the  assertion  that  X  <  Y  <  X . 
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equations  with  increasing  degree  of  degeneracy.  In  particular,  consider  the  two-dimensional 
equation 

pt  +  (9X1  +  dX2)A(p)  —  (dxi  —  d X2)~B(p)  =  0. 

If  we  set  A  =  0,  the  equation  has  a  strong,  rank-one  parabolic  degeneracy  with  no  regularizing 
effect  coming  from  its  purely  diffusion  part,  since  (b(u)£,£)  =  0,  V£i  —  £2  =  0,  indicating 
the  persistence  of  steady  oscillations,  po(x  +  y );  moreover,  if  B  —  0  then  the  equation  has  no 
regularization  coming  from  its  purely  convection  part,  since  a(u)  •  £  =  0,  V£i  +£2  =  0  indicates 
the  persistence  of  steady  oscillations  p0(x  —  y).  Nevertheless,  the  combined  convection- diffusion 
with  A(p)  ~  pl+1  and  B(p )  ~  \p\np  does  have  IW1 -regularizing  effect  of  order  s  <  6/ (2  +  2 n—£) 
for  n  >  2d,  consult  corollary  4.5  below. 

Finally,  in  Section  5  we  consider  degenerate  elliptic  equations,  —J2df.XkBjk(p)  =  S(p).  As¬ 
suming  that  the  non-degeneracy  condition  (1.4)  holds,  then  the  kinetic  formulation  of  bounded 
solutions  for  such  equations,  p  G  IVs'1  (D),  have  interior  regularity  of  order  s  <  min(o;,  2a/{2a  + 
1)).  We  conclude  by  noting  that  it  is  possible  to  adapt  our  arguments  in  more  general  setups, 
for  example,  when  a  suitable  source  term  S(p)  it  added  to  the  time  dependent  problem,  when 
dealing  with  degenerate  temporal  fluxes,  dtA0(p),  A'0  >  0,  or  when  lower-order  convective  terms, 
V.x  •  A(p),  are  added  in  the  elliptic  case. 

2.  The  averaging  lemma 

We  are  concerned  with  the  regularity  of  averages  of  solutions  for  differential  equations  of  the 
form 

(2.1)  C(V„v)f  =  A+V  A*:=(— A +A 

Here,  /  =  f(x,  v)  G  WffffWf  x  !+)  and  g  =  g(x,  v )  G  Lqloc(Rx  x  M„)  are  real-valued  functions  of 
the  spatial  variables  x  =  (aq, . . .  ,  xf)  G  and  an  additional  parameter  dgI,  called  velocity 
by  analogy  with  the  kinetic  framework,  and  jC(Vx,v)  is  a  differential  operator  on  of  order 
<  k,  whose  coefficients  are  smooth  functions  of  v. 

The  velocity  averaging  lemma  asserts  that  if  £(-,v)  is  nondegenerate  in  the  sense  that  its 
null  set  is  sufficiently  small  —  to  be  made  precise  below,  then  the  u-moments  of  f(x,  •), 

7{x)  ■=  I  f{x,v)</>(v)dv,  0GC“, 

J  V 

are  smoother  than  the  usual  regularity  associated  with  the  data  of  f(x,-)  and  g(x,-).  That 
is,  by  averaging  over  the  so-called  microscopic  u-variable,  there  is  a  gain  of  regularity  in  the 
macroscopic  x- variables.  There  is  a  relatively  short  yet  intense  history  of  such  results,  motivated 
by  kinetic  models  such  as  Boltzmann,  Vlasov,  radiative  transfer  and  similar  equations  where 
the  v- moments  of  /  represent  macroscopic  quantities  of  interest.  We  refer  to  the  early  works  of 
Agoshkov  [Ag84]  and  Golse,  Lions,  Perthame  and  Sentis,  [GPS85,  GLPS88]  treating  first-order 
transport  operators  with  /,  g  integrability  of  order  p  =  q  >  1.  The  work  of  DiPerna,  Lions  and 
Meyer,  [DLM91],  provided  the  first  treatment  of  the  general  case  p  f  q,  followed  by  Bezard, 
[Be94] ,  their  optimality  in  [Li95]  and  an  optimal  Besov  regularity  result  of  DeVore  and  Petrova, 
[DVP00]. 

Extensions  to  more  general  streaming  operators  were  treated  by  DiPerna  and  Lions  in 
[DL89a,  DL89b]  with  applications  to  Boltzmann  and  Vlasov-Maxwcll  equations,  and  by  Li¬ 
ons,  Perthame  and  Tadmor  in  [LPT94a]  with  applications  to  nonlinear  conservation  laws 
and  related  parabolic  equations.  Gerard,  [Ge90]  together  with  Golse  [GG92]  provided  an 
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L2-treatment  of  general  differential  operators.  A  different  line  of  extensions  consists  of  ve¬ 
locity  lemmata  which  take  into  account  different  orders  of  integrability  in  x  and  v,  leading  to 
sharper  velocity  regularization  results  sought  in  various  application.  We  refer  to  the  results 
in  [DLM91]  for  /  G  Lpl(Rv,  LP2(M^)),  g  G  Lqi  (Rv ,  Lq2  (R^))  and  to  Jabin  and  Vega  [JV03]  for 
/  G  WNl’pi(RVl  LP2(R^)),  g  G  L92(Rf)).  Westdickenberg,  [We02],  analyzed  a  general 

case  of  the  form  /  G  B^^fRf ,  Lri(Rv)),  g  G  Bf  <n (R^ ,  L7"2 (R„ ) ) .  Jabin,  Perthame  and  Vega, 
[JP02],  [JV04]  used  a  mixed  integrability  of  /  G  LP(R^,  kPAri,p(M.l,)),  g  G  Lq(R^.,WN2,q(Rv))  to 
improve  the  regularizing  results  for  nonlinear  conservation  laws  [LPT94a] ,  [LPT94b] ,  Ginzburg- 
Landau,  and  other  nonlinear  models.  Golse  and  Saint-Raymond,  [GSR02]  have  shown  that  a 
minimal  requirement  of  equi-integrability,  say  /  G  L^R^,  L®(RV)),  g  G  L^R^  x  R„)  measured 
in  Orlicz  space  L #  with  super-linear  <F,  is  sufficient  for  relative  compactness  of  the  averages  /, 
which  otherwise  might  fail  for  mere  L1 -integrability,  [GLPS88]. 

The  derivation  of  velocity  averaging  in  the  above  works  was  accomplished  by  various  methods. 
The  main  approach,  which  we  use  below,  is  based  on  decomposition  in  Fourier  space,  carefully 
tracking  f(£,v)  in  the  “elliptic”  region  where  [v  :  £(i£,v)  ^  0},  and  the  complement  region 
which  is  made  sufficiently  small  by  a  non-degeneracy  assumption.  Other  approaches  include  the 
use  microlocal  defect  measures  and  H-measures  [Ge90,  Ta90],  wavelet  decomposition  [DVP00] 
and  “real-space  methods”  —  in  time  [Va99,  BD99],  and  in  space-time  using  Radon  transform 
[ChOO,  We02],  X-transform  [JV03,  JV04]  and  duality-based  dispersion  estimates  [GSR02],  Al¬ 
most  all  of  these  results  are  devoted  to  the  phenomena  of  velocity  averaging  in  the  context  of 
transport  equations,  k  —  1, 

Our  study  of  velocity  averaging  applies  to  a  large  class  of  £’s,  satisfying  the  so-called  trun¬ 
cation  property,  in  Section  2.4  we  show  that  all  £’s  of  order  k  <  2  satisfy  this  truncation 
property.  In  particular,  we  improve  the  regularity  statement  for  first-order  velocity  averaging 
and  extend  the  various  velocity  averaging  results  of  the  works  above  from  first-order  transport 
to  general  second-order  transport-diffusion  and  elliptic  equations.  The  results  are  summarized 
in  the  averaging  lemmata  2.1  and  2.2  for  homogeneous  symbols  and  in  averaging  lemma  2.3 
for  general,  truncation-property-satisfying  C s.  Our  derivation  is  carried  out  in  Fourier  space 
using  Littlewood-Paley  decompositions  of  /  G  x  R„)  and  g  G  Lloc(Rf.  x  R„).  To  avoid 

an  overload  of  indices,  we  leave  for  future  work  possible  extensions  for  more  general  data  with 
mixed  (x,  v)-integrability  of  /  and  g. 

2.1.  The  truncation  property.  We  now  come  to  a  fundamental  definition. 

Definition  2.1.  Let  rn(f  )  be  a  complex- valued  Fourier  multiplier.  We  say  that  m  has  the  trun¬ 
cation  property  if,  for  any  locally-supported  bump  function  if  on  C  and  any  1  <  p  <  oo, 
the  multiplier  with  symbol  ?/?(m(£)/d)  is  an  lA-multiplier  uniformly  in  5  >  0,  that  is,  its  Lp- 
multiplier  norm  depends  solely  on  the  support  and  size  of  if  (for  some  large  £  which  may 
depend  on  m)  but  otherwise  is  independent  of  5. 

In  Section  2.4  we  will  describe  some  examples  of  multipliers  with  this  truncation  property. 
Equipped  with  this  notion  of  a  truncation  property,  we  turn  to  discuss  the  Lp-size  of  parame¬ 
terized  multipliers.  Let 

M^f(x,  v)  :=  Tfflif(m(y^  —  )/(£w),  where  f(£,v)  =  Fxf{£,v)  :=  f 

v  o  J  J 


e~2^  xf(x,  v )  dx 
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denote  the  operator  associated  with  the  complex-valued  multiplier  truncated  at  level 

6.  Our  derivation  of  the  averaging  lemmata  below  is  based  on  the  following  straightforward 
estimate  for  such  parameterized  multipliers. 

Lemma  2.1  (Basic  Lp  estimate).  Let  /  be  a  finite  interval,  I  C  and  assume  m(£,v)  satisfies 
the  truncation  property  uniformly  in  v  G  I.  Let  1  <  p  <  2.  Let  M.^  denote  the  velocity- 
averaged  Fourier  multiplier 

M^f(x)  :  =  j M^f(x,v)  dv  =  J J7~1ip(m(^,v)/S)Jrxf(x,v)  dv. 

For  each  (etf  and  <5  >  0  let  5)  C  /  be  the  velocity  set 

nm({;  5)  :=  {v  e  I  :  V'>  e  supp  V^j. 

Then  we  have  the  Lp  multiplier  estimate 

(2-2)  ll-Ady/(a;)||Lp(Rd)  <  sup  <J)|^  •  \\f\\Lp(RdxRv) 

Proof.  For  p  —  2,  the  claim  (2.2)  follows  from  Plancherel’s  theorem  and  Cauchy-Schwarz  in¬ 
equality,  while  for  p  close  to  1+  the  claim  follows  (ignoring  the  bounded  |Qm(£;  5)  | -factor)  from 
the  assumption  that  m(£,v)  satishes  the  truncation  property  uniformly  in  v  (and  in  fact,  the 
end  point  p  —  1,  with  the  usual  Lt 1  -replacement  of  L1,  can  be  treated  by  a  refined  argument 
along  the  lines  of  [DLM91]).  The  general  case  of  1  <  p  <  2  follows  by  interpolation.  □ 

Remark  2.1.  Clearly,  if  m(£,  •)  and  c(£)  are  //-multipliers  then  so  is  their  product,  and  in 
particular,  if  m  has  the  truncation  property  then  (2.2)  applies  for  'if>(m(f,v)/S)c(f). 

2.2.  Averaging  lemma  for  homogeneous  symbols.  We  will  present  several  versions  of  the 
averaging  lemma.  The  later  versions  will  supercede  the  former,  but  for  pedagogical  reasons 
we  will  start  with  the  simpler  case  of  homogeneous  symbols.  In  this  case  it  is  convenient  to 
use  polar  co-ordinates  f  =  |£|£',  where  f  G  S^-1  is  defined  for  all  non-zero  frequencies  f  by 

£':=£/l£l- 

We  begin  with 

Averaging  Lemma  2.1.  Let  1  <  q  <  2  and  let  g  G  Lgloc(Rd  x  M„)  if  q  >  1  or  let  g  be  a  locally 
bounded  measure,  g  G  -M(Mf  x  M„)  if  q  =  1.  Let  rp  N  >  0  arid  let  f  G  WffiffMf.  x  Rv),  a  > 
0,  1  <  p  <  2  solves  the  equation 

(2.3)  C{Vx,v)f{x,v)  =  A l^g{x,v)  in  V'(Rdx  x  Rv). 

Here,  £(VX,-)  is  a  differential  operator  with  sufficiently  smooth  coefficients,  £(-,v)  G  Cn,£>0, 
and  let  £(i£,v)  be  the  corresponding  symbol.  We  assume  that  £(if,v)  is  homogeneous  in  f  of 
order  k,k  >  a  +  rj,  that  the  modified  symbol  £{if  ,v)  obeys  the  truncation  property  uniformly 
in  v,  and  that  it  is  nondegenerate  in  the  sense  that  there  exists  an  a,  0  <  a  <  (N  +  l)q' ,  such 
that 

sup  |n£(£;<5)|  <  5a,  f lc(f;5)  :=  {v  G  /  :  |£(if,u)|  <  <*}. 

$eRd:|5|  =  l  1  J 


(2.4) 
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Then,  there  exist  9  =  9a  G  (0, 1)  and  sa  =  s(9a )  >  a  such  that  for  all  bump  functions 

locV^xJ  f°r 


4>  G  Cq°(I),  the  averages  f(x)  :=  f  f(x,v)f>(v)dv  belong  to  the  Sobolev  space 


all  s  G  (a,  sa)  arid  the  following  estimate  holds 

(2-5)  H/||w,s0£(Rg)  ~  \\d  II  xR„)  '  Wf\\wf0f(RixRv] 
Here  6  =  9a(p,  q,  N )  and  r  are  given  by 

a/p'  1 


s  G  (a,  sa),  sa  :=  (1  -  9a)a  +  9a(k  -  rj). 


(2.6) 


9  = 


a(l/p' —  1/q')  +  N  +  l1  r 


1-9  9 

- 1- 

p  q 


o  <  9  <  1. 


Remark  2.2.  It  would  be  more  natural  to  assume  that  the  symbol  C(if,v)  itself,  rather  than 
the  modified  symbol  £(i£',v),  obeyed  the  truncation  property,  as  it  is  typically  easier  to  verify 
the  truncation  property  for  the  unmodified  symbol.  Indeed  when  we  turn  to  more  advanced 
versions  of  the  averaging  lemma  (which  rely  on  Littlcwood-Paley  theory,  and  which  do  not 
assume  homogeneity)  we  will  work  with  the  truncation  property  for  the  unmodified  symbol. 
However  we  choose  to  work  here  with  the  modified  symbol  as  it  simplifies  the  argument  slightly. 


Proof.  We  start  with  a  smooth  partition  of  unity,  1  =  )T)  'ifj(2~iz)  such  that  ifo  is  a  bump 
function  supported  inside  the  disc  \z\  <  2  and  the  other  'ipj’s  are  bump  functions  supported  on 
the  annulus  1/2  <  \z\  <  2  (we  note  in  passing  that  the  other  yl/s  can  be  taken  to  be  equal,  so 
the  index  j  merely  serves  to  signal  their  ‘action’  on  the  shells,  2J_1  <  \z\  <  2J,+1).  We  set 


fj(x,v)  :=  J  =  0,1,2...  , 

recalling  that  f  :=  £/|£|,  and  we  consider  the  corresponding  decomposition  f  —  fo  +  Ylj>i  fj- 
We  distinguish  between  two  pieces,  /  =  f^+f^  where  /(°)  :=  /0  and  /d)  :=  Y2j>i  fj-  Observe 

that  the  w-support  of  /(°)  is  restricted  to  the  degenerate  set  f/c(£;t>)  whereas  /d)  =  Ylj>ifj 
offers  a  decomposition  of  the  non-degenerate  complement,  £!£(£;  5).  The  free  parameter  6  is  to 
be  chosen  later. 

We  start  by  noting  that  /(0)  =  /0  is  associated  with  the  multiplier  rf0(£(if',v)/5).  Since 
£(if',v)  satisfies  the  truncation  property,  we  can  use  Lemma  2.1  and  non-degeneracy  assump¬ 
tion  (2.4),  to  obtain 


(2.7) 


ll/(0)ll 


wpp 

loc 


sup 

1*1=1 


1  tv’ 


(Rg  xR„) 


xR„)- 


We  turn  to  the  other  averages,  fj,  j  >  1  which  make  /d).  Since  £(if,-)  is  homogeneous  of 
order  k,  equation  (2.3)  states  that 


lei 


M-'m.v)  =  A 


and  thus, 


a  rvi  =  x-7^-1  / 1 


x  2i5  x  ;  V  2 i8 


<%9(Z,v)<i>(v)dv,  j  =  1,2,..., 


where  'ifj(z)  :=  ifj(z)/z  is  a  bump  function  much  like  ipj  is.  Integration  by  parts  then  yields 


(2.8)  A ZT’fj  = 


i5)N+1 


(2i5) 


T, 


-i 


iN)  f£(if',v) 


J  V  2  35 


£v(f&  v)  </(£•  v)o(i.')dv  - 


+  lower  order  or  similar  terms. 
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We  can  safely  neglect  the  lower  or  similar  order  terms  which  involve  (powers  of)  the  bounded 
multipliers  ,  v) ,  l  <  N  and  we  focus  on  the  leading  term  in  (2.8),  associated  with  the 

multipliers 


(2.9) 


-{N) 


(£(%',  y) 

V  2  36 


By  Hormander-Mikhlin  or  Marcinkeiwicz  multiplier  theorems,  C™  are  bounded  multipli¬ 

ers  and  hence  (2.8)  are  upper-bounded  by 


IIA^V, 


I  LI, 


(2  SS) 


N+ 1 


Llo, 


ffere  A4j  =  A4~(n)  are  the  Fourier  multipliers  with  symbol  i/jj  (£(i£',  v)/2q5). 

ifjj 

We  now  fix  q  >  1.  By  our  assumption,  £(i£',  v)  satisfies  the  truncation  property,  and  Lemma 
2.1  implies  (Rd\  ;$  (2q+15)a^ql \\g\\Lq  (R<q-  Adding  together  all  the  f/s,  we  hnd  that 

/(1)  =  i  fj  satisfies 


(2.10)  IIATV^IIi*  «*) 


3>  1 


(2^) 


N+l 


(2 J+1h)- 


ifocC 


5i‘ 


—  (JV+l) 


Ll„c( 


Thus,  if  we  fix  t  >  0  and  choose  5  to  equilibrate  the  bounds  in  (2.7)  and  (2.10), 


s«W-W+»+  1  ~  t\\g(x,v)\\L,J\\nx,v) 


I  "ST* 


then  this  tells  us  that 


_ inf 

/(°>+/(  1)=7  L 


Wf{0)Ww?*m)  +  t\\f{1)\\wlk-^(Ri)  ~  f  '  lb(xw)||^oc(RgxR,)  '  11/(^011 


1—9 

W^(RZxRvy 


and  the  desired  llq^-bound  follows  for  s  <  (1  —  9) a  +  9{k  —  rj)  with  6  given  in  (2.6).  The 
remaining  case  of  q  =  1  can  be  converted  into  the  previous  situation  using  Sobolev  embedding. 
In  this  case,  g  being  a  measure,  it  belongs  to  W~e,qe  for  all  (e,  qe)  such  that 


g  G  W~e’q% 


d  +  2 

Qe 


<  e  <1  <  qe  < 


d  +  2 
d  +  1’ 


and  hence  (2.5)  applies  for  s  <  (l—9)a  +  9[k  —  g  —  e)  and  6  =  9a(p,  qe,  N);  we  then  let  e  approach 
0+  so  that  qe  approaches  arbitrarily  close  to  1+  to  recover  (2.5)  with  sa  and  9a(p,  1,77).  □ 


Remark  2.3.  As  an  example  consider  a  (possibly,  pseudo-) differential  operator  £(V x,  •)  of  order 
k  and  let  f(x,v)  G  such  that  £(Vx,v)f  G  W°Q~k+1’2.  Assume  that  Cv(i£,v)  ^  0  so  that 
the  nondegeneracy  condition  (2.4)  holds  with  a  —  1.  Application  of  the  averaging  lemma  2.1 
with  p  =  q  =  2,N  =  0  and  rj  —  k  —  a  —  1  then  yields  the  gain  of  half  a  derivative,  f{x)  G 
with  s  <  cr/2  +  (cr  +  l)/2  =  a  +  1/2,  in  agreement  with  [GG92,  Theorem  2.1] .  The  main  aspect 
here  is  going  beyond  the  L2-framework,  while  allowing  for  general  and  possibly  different  orders 
of  integrability,  (f,g)  G  (Wa'p,  Lq). 


Remark  2.4.  The  limiting  case  of  the  interpolation  estimate  (2.5),  9  —  1,  s  —  sa,  corresponds 
to  Besov  regularity  /  G  This  regularity  can  be  worked  out  using  a  more  precise 

bookkeeping  of  the  Littlewood-Paley  blocks.  For  the  transport  case,  k  —  1,  it  was  carried 
out  first  in  [DLM91,  Theorem  3],  improved  in  [Be94]  and  a  final  refinement  with  a  secondary 
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index  t  =  p  can  be  found  in  [DVPOO].  This  limiting  case  is  encountered  in  the  particular 
situation  when  k  =  77,  so  that  the  interval  (a,  sa )  ‘survives’  at  9=1.  Here,  one  cannot 
expect  for  regularizing  effect,  but  there  is  a  persistence  of  relative  compactness  of  the  mapping 
g(x,v)  ^  J(x ),  [PS98], 


To  gain  a  better  insight  into  the  last  averaging  lemma,  we  focus  our  attention  on  the  case 
where  /  is  a  Wfififf- -solution  of 

(2.11)  £(' V.r,  v)f(x,  v)  =  dvg(x,  v),  g  G  Lqloc(R*  x  Rv),  /  G  Wfifc(Rdx  x  Rv), 


corresponding  to  the  special  case  77  =  0  and  N  =  1  in  (2.3).  This  case  will  suffice  to  cover  all 
the  single- valued  applications  we  have  in  mind  for  the  discussion  in  Sections  3,  4  and  5  without 
the  burden  of  carrying  out  an  excessive  amount  of  indices.  The  averaging  lemma  (2.1)  implies 
that  f(x)  has  Sobolev  regularity  of  order  s  <  (1  —  9)a  +  9k, 


(2.12) 


LL 


«) 


i-e 

lp1oc 


6  =  6n  := 


a/p' 


a^/p1-  W)  +  2' 


The  last  regularity  statement  can  be  improved.  To  this  end,  we  revisit  the  dyadic  multipliers 


. — -IN) 

in  (2.9),  7 />j  v)/2:’S)£,(f (i£',v).  The  key  observation  is  that  Cv(i£f  ,v)  acts  only  on  the 

subset  of  u’s  —  those  which  belong  to  v  G  2-?+1<5).  Linking  the  size  of  v)  to  that 

of  C(i£/,v),  we  arrive  at  the  following  improved  averaging  regularity  lemma. 


Averaging  Lemma  2.2.  Let  f  G  W££(Rd  x  R„),  o  >  0, 1  <  p  <  2,  solves  the  equation 

,91^  rm  \f(  \  a  (  \  1  <  g  <  2, 

(2.13)  C(Vx,v)f(x,v)  =  dvg(x,v),  g  G  j  x  q  =  L 


Let  £(*£,  v )  be  the  corresponding  symbol.  We  assume  that  C{i £,  v)  is  homogeneous  in  i £  of  order 
k,  k  >  a,  that  the  modified  symbol  C[if  ,v)  satisfies  the  truncation  property  uniformly  in  v,  and 
that  it  is  nondegenerate  in  the  sense  that  there  exists  an  a,  0  <  a  <  q' ,  such  that  (2.Jfi)  holds. 
Moreover,  assume  that 


(2.14)  3/7  G  [0, 1]  s.t.  sup  sup  | £„(*£, v)\  <  £|c(£;<5) 

kl=i  ven£(W) 


{v  e  I  :  |£(if,u)|  <<S}. 


Then,  there  exist  6 
(2.15) 


9a  G  (0, 1)  and  r  given  by 

Q  _ _ 

a(l/p'  :-  l/q')  +  2  -  p 


1 

r 


1-9 

p 


such  that  for  all  bump  functions  G  C/fi(L),  the  averages  f(x)  :=  f  f(x,v)(f>(v)dv  belong  to  the 
Sobolev  space  Wls^(R/c)  for  all  s  G  (a,  sa),  sa  =  (1  —  9a)cr  +  9ak  and  (2.5)  holds. 


For  the  proof  we  revisit  (2.8)  with  77  =  0,  N  =  1, 

Axfj  =  d^3  (C^2js  ^  )  Cv v)(f>(v)dv . 

Its  bound  in  (2.10)  can  now  be  improved  by  the  extra  factor  of  (2 which  follows  from 
(2.14),  yielding 

\\A-xf j\\l1oc(r£)  ~  (2iS) 2 

j>  1  v  ' 

and  we  conclude  by  arguing  along  the  lines  of  the  averaging  lemma  2.1.  □ 


(2 j+15) 


■b+M 

q 


LL 


LL 


«)1 


10 


EITAN  TADMOR  AND  TERENCE  TAO 


Remark  2.5.  In  the  generic  case  of  a  symbol  C(i£,v)  which  is  analytic  in  if  uniformly  in  v,  the 
‘degeneracy  of  order  a’  in  (2.4)  implies  that  sup  {|£„(i£,  v)|  :  |£|  =  1,  |£(i£,u)|  <  5}  <  h1_". 
Thus,  (2.14)  holds  with  p  —  1  —  a  and  in  this  case,  the  velocity  averaging  holds  with  Sobolev- 
regularity  exponent 


f(x)  e 


s  <  sa  :=  (1  —  0a)cr  +  6ak,  9a  := 


<*/p' 

a(l/p'  +  l/q)  +  T 


Remark  2.6.  The  statement  of  the  averaging  lemma  2.2  can  be  extended  for  /  G  LP(R^  x  R„) 
and  g  G  L^oc(M dx  x  R„)  in  the  full  range  of  1  <  p  <  oo  and  1  <  q  <  oo.  The  dual  claim  to  (2.2), 
based  on  (L2,  BMO )  interpolation  for  2  <  p  <  oo  reads 


(2.16) 


\\M^f(x)\\LP{^d)  <  sup  |fim(f;  S)\p  •  || /|| Lp(s.d xR„)i  2  <  p  <  oo. 
C 


This  yields  the  same  IT  ^-regularity  as  before 


(2.17) 


LL 


)  '  ll/ll  VK,‘T’p(RgxR„)’  S  ^  (T  ®a)a  +  @otk) 


with  6  =  0a(p,  q)  given  by 
a/p 


^  ^  _  Q  Q 

0  :=  -777= — 7 '/-\  ,  o - ’  _  :=  +  =’  P  '■=  ma x(p,p')>  9  :=  max(g,  g') 

a(l/p  —  1/g)  +  2  —  p  r  p  q 


2.3.  An  averaging  lemma  for  general  symbols.  We  now  turn  our  attention  to  averages 
involving  general,  not  necessarily  homogeneous  symbols;  as  such  the  polar  co-ordinate  repre¬ 
sentation  f  =  is  no  longer  useful  and  will  be  discarded.  We  focus  on  equations  of  the  form 
£(yx,v)f(x,v)  =  dvg(x,v),  corresponding  to  (rj ,N)  =  (0,1)  in  (2.3). 


Averaging  Lemma  2.3.  Let  f  G  IT/7 (®7  x  R„),  a  >  0, 1  <  p  <  2,  solves  the  equation 


(2.18) 


£(Vx,v)f(x,v)  =  d Vg(x,v), 


f  Li(Ri  x  R„),  1  <  q  <  2, 

^(^xt0),  g=l. 


Let  £(if,  •)  be  the  corresponding  symbol  of  degree  <  k  with  sufficiently  smooth  v-dependent 
coefficients  and  assume  it  obeys  the  truncation  property.  Denote 

uc{J]6):=  sup  |Uc(£;<f)|,  8)  :=  {v  :  |£(i£,  v)|  <  <5} 

£eRd:| 

and  suppose  the  following  non- degeneracy  condition  holds 

(2.19)  3o,/3>0  s.t.  j  ,  Vh  >  0,  J  >  1. 

Moreover,  assume  that 

(2.20)  3A  >  0  and  [i  G  [0,1]  s.t.  sup  sup  |£„(i£,  v)|  <  J;3AT. 

\t\~J  venc(t;S) 


Then,  for  all  bump  functions  /  G  C™(I),  the  average  f(x)  :=  f  f(x,v)<f(v)dv  belongs  to  the 
Sobolev  space  W^( R^)  for  s  G  (cr,saj 3)  and  the  following  estimate  holds 


(2.21) 


f(x,v)(f>(v)  dv 


wfM) 


Wf(x,v)\\w^(RixRv)  +  \\9{xD’)\\Lloc(RdxRvy 
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Here,  sa^  :=  (1  —  9a)a  +  0a/3{2  —  p  —  A)  where  9  =  9a  and  r  are  given  by 


(2.22) 


a/p'  1  1  —  9  9 

a(l/p’  —  1/q')  +  2  —  /i’  r  p  q ’ 


0  <  0  <  1. 


Remark  2.7.  How  does  the  last  averaging  lemma  compare  with  the  previous  ones?  We  note 
that  since  deg£y(i£,  •)  <  k,  then  the  additional  assumption  (2.20)  always  holds  with  p  =  0  and 
f3 A  =  k.  Hence,  the  averaging  lemma  2.3  with  just  the  non-degeneracy  condition  (2.19)  yields 
the  regularity  f(x)  G  W'^(R^)  of  (the  reduced)  order  s  G  (<j,satg),  where 


(2.23) 


sa,p  ■=  (1  -  9a)a  +  9a(2f3  -  k ), 


a/p' 

a(l/p’-  1/q') +  2’ 


Now,  if  in  particular,  £{i £,  •)  is  homogeneous  of  degree  k ,  then 

u  £■(&'■>  d)  =  ^  £  -=  jlj! 


this  shows  that  if  the  non-degeneracy  condition  (2.4)  of  averaging  lemma  2.1  holds,  d)  ~ 

(5/Jfc)°,  then  it  implies  (2.19)  with  (3  =  k,  and  we  recover  the  homogeneous  averaging  lemma 
2.1  with  (77,  N)  =  (0, 1),  namely,  the  averages  /  gain  regularity  of  order  s  <  {l—9)a+9{2j3—k)  = 
(1  —  9)cr  +  9k.  The  only  difference  is  that  now  the  truncation  property  is  assumed  on  the 
unmodified  symbol  £(i£,u)  rather  than  the  modified  one  C(i^',v). 

As  for  the  averaging  lemma  2.2,  we  first  note  that  in  the  generic  case  of  an  homogeneous 
£(-,v),  the  additional  assumption  (2.20)  holds  with  A  =  a  and  p  =  1  —  a, 

(2.24)  sup  sup  \£v(if/,v)\  < 

|£|~J  {vel:  |£(«, «)!<$} 


Indeed,  all  the  homogeneous  examples  discussed  in  S3ections  3,4  and  5  below,  employ  the 
averaging  lemma  2.2  with  these  parameters  which  yield  IT '^-regularity  of  order  s  <  (1  — 
9a)cr  +  9a(3, 


(2.25) 


/(*)  e  WC(Rj),  s  <  (1  -  ea)a  +  M,  0a 


a/p' 

a(l  jp'  +  1/q)  +  1 


In  the  particular  case  of  £  being  homogeneous  of  order  k  then  /3  =  k  and  we  recover  the 
averaging  lemma  2.2  (except  that  the  truncation  hypothesis  is  now  assumed  on  the  unmodified 
symbol) . 


Proof.  We  begin  by  noting  that  we  can  safely  replace  the  non-degeneracy  condition  (2.19)  with 
a  slightly  weaker  one,  namely 

/(51+E\a 

(2.26)  3a,  P  >  0,  and  e  >  0  s.t.  5)  ~  \  JJTJ  ’  V5  >  0,  J  >  1, 

and  still  retain  the  same  gained  of  regularity  of  order  s  <  sa This  can  be  achieved  by 
replacing  the  values  of  a  in  (2.19)  by  a/(l  +  e)  and  then  absorbing  £  into  a  slightly  smaller 
order  of  regularity,  (1  —  e)<sQj(a  dictated  by  9ap i+£).  The  extra  e-power  of  5  will  be  needed  below 
to  insure  simple  summability,  which  probably  could  be  eliminated  by  a  more  refined  argument 
involving  Besov  spaces,  along  the  lines  of  [DLM91]. 

Next,  we  break  up  /  into  Littlewood-Paley  pieces, 

f  —  fo+  ^  fj, 

dyadic  1 
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so  that  fj(£,v),  the  spatial  Fourier  transform  of  fj(x,v),  is  supported  for  frequencies  |£|  ~  J, 
and  /o  has  support  in  |£|  <  1.  Since  f0  is  a  smooth  average  of  /  at  unit  scales,  the  contribution 
of  /o  is  easily  seen  to  be  acceptable.  By  giving  up  an  £  in  the  index  sajp  we  may  thus  reduce 
(2.21)  to  a  single  value  of  J.  It  thus  suffices  to  show  that 


JS  ||/j||l[  (Rd)  ~  Wfj  II  Tfoc(Ri  xl„)  + 


£?oe(*i> 


„)>  fj  '■=  /  fj(x,v)(f)(v)dv, 


for  each  J  >  1. 

Fix  J.  Because  of  the  local  nature  of  Littlewood-Paley  projections  when  J  >  1  we  may 
replace  the  localized  Lp  norms  with  global  norms.  Actually  we  may  replace  Lr  by  weak  Lr 
since  we  may  pay  another  £  in  the  index  saip  to  improve  this.  By  the  duality  of  weak  Lr  and 
Lr>'1  it  thus  suffices  to  show  that 


(2.27) 


\(fj,XE)\  &  J-\E\1'r' 


for  all  sets  E  C  of  finite  measure,  where  we  have  normalized  ||/j||lp(r^xr„)  ~  J  a  and 


Llc(RixRv) 


<  1. 


We  now  decompose  the  action  in  u-space  of  each  of  the  Littlewood-Paley  pieces  (rather  than  a 
decomposition  /  itself  used  in  lemma  2.1), 


«*,»)=  £ 

dyadic  5's~Jk 


'£(Vx,v) 


'£{i£,v) 


J-r  . 


Here,  'ijj(z)  is  a  bump  function  on  C  supported  on  the  region  \z\  ~  1.  It  suffices  to  estimate 

'C(Wx,v) 


(2.28) 


0 


fj{~i  v)4>{v)dv,  xe 


with  a  summable  decay  as  the  dyadic  <5  — >  0,  so  that  (2.27)  holds. 

By  our  assumption,  £(i£,  v)  satisfies  the  truncation  property  uniformly  in  v,  hence  by  Lemma 
2.1  we  see  that  for  all  1  <  p  <  2, 


(2.29) 


0 


5 


fj(x,  v)4>{y)dv 


Uc{J’i  d)l^P  II  f J  II  Tf  (R^  xR„)  • 


From  (2.29)  and  Holder  we  may  thus  estimate  (2.28)  by 
(2.30) 


0(— ^f^)fj(-,V)(l>(V)dv,xE)  ZJ  auc(J-,5)1/pl\E\1/p' . 


On  the  other  hand,  thanks  to  equation  (2.18)  we  can  write 


0 


c(yx,v) 

5 


fj(x,v )  =  0 


£(V*,u)\l  d 


5 


S  dv 


9j(xi  v) 


where  ip(z)  :=  ip(z)/z  and  the  gf  s  are  the  corresponding  Littlewood-Paley  dyadic  pieces  of  g. 
We  thus  have 


(2.31) 


0 


£(Vx,f) 

5 


fj(x,  v)(f>(v)dv  =  -  /  0 


C{Vx,v)\  d 


5 


dv 


9j(x,  v)4>(v)dv. 
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We  now  integrate  by  parts  to  move  the  d/dv  derivative  somewhere  else.  We  will  assume  that 
the  derivative  hits  i()(£(S7x,  v)/§),  as  the  case  when  the  derivative  hits  the  bump  function  (f>(v) 
is  much  better.  We  are  thus  led  to  estimate 

(2-32)  h  (f  ~ ~ )  Cv^x'  v)9j(xi  v)<j>(v)dv,  Xe)  ■ 

Since  gj  is  localized  to  frequencies  ~  J,  then  by  (2.20),  the  multiplier  Cv(i((,v)  acts  like  a 
constant  of  order  O  (J ,3X  5^) .  Also,  i(z  is  a  bump  function  much  like  i/j.  Thus  we  may  modify 
(2.29)  —  with  p  replaced  by  q  (assuming  that  q  >  1  and  using  the  modified  argument  for  the 
case  q  =  1  as  before),  if}  replaced  by  ^z,  and  fj  replaced  by  £v{Vx,v)gj,  to  estimate  (2.31)  by 

(2.33)  |(  J  il>(C{W*,V))fj(-,v)(l)(v)dv,XE)\  Z  r^J^V(^5)W|A|W. 

Interpolating  this  bound  with  (2.30),  we  may  bound  (2.28)  by 

<$-*(2-**)  ( J;  5)  W'  I E  |  l'r' . 

The  parameterization  in  (2.22)  dictates  a/r'  =  9(2  —  p).  Finally,  we  put  the  extra  £  power  into 
use:  by  (2.26),  uc(J]  S)1^'  £  (£1+e ll>  and  hence  the  last  quantity  is  bounded  by 
^(2-^j-^i^i/p  Sq/3  =  (!  _  Q)a  -  9(3(2  -  g  -  A). 

Summing  in  <5  and  using  the  hypothesis  that  s  <  sa^  we  obtain  (2.27).  □ 

2.4.  Velocity  averaging  for  first  and  second  order  symbols.  To  apply  the  velocity  aver¬ 
aging  2.3  we  need  to  find  out  which  multipliers  m(£)  have  the  truncation  property.  Fortunately, 
there  are  a  large  classes  of  such  multipliers. 

First  of  all,  it  is  clear  that  the  multipliers  m(£)  =  £  •  e\  and  m(£)  =  |^|2  have  the  truncation 
property,  as  in  these  cases  the  Fourier  multipliers  are  just  convolutions  with  finite  measures. 
Now,  observe  that  if  m(£)  has  the  truncation  property,  then  so  does  m(L( £))  for  any  invertible 
linear  transformation  L  on  with  a  bound  which  is  uniform  in  L.  This  is  because  the  Lp 
multiplier  class  is  invariant  under  linear  transformations. 

Because  of  this,  we  see  that  the  multipliers 

™>i(0  =  a(u)  -i£ 

and 

^2(0  =  (b(u)f,f) 

have  the  truncation  property  uniformly  in  v,  where  a(u)  are  arbitrary  real  coefficients,  and  b(u) 
is  an  arbitrary  elliptic  bilinear  form  with  real  coefficients. 

From  the  Hormander-Mikhlin  or  Marcinkeiwicz  multiplier  theorems  and  the  linear  transfor¬ 
mation  argument  one  can  also  show  that  mi((')  has  the  truncation  property  uniformly  in  v. 
These  arguments  go  back  to  the  discussion  of  [DLM91].  The  situation  with  n^^)  is  less  clear, 
but  fortunately  we  will  not  need  to  verify  that  these  second  order  modified  multipliers  obey 
the  truncation  property  since  our  averaging  lemmata  also  work  with  a  truncation  property 
hypothesis  on  the  unmodified  multiplier. 

Now  we  observe  that  if  mi(£),  m2( £)  are  real  multipliers  with  the  truncation  property,  then 
the  complex  multiplier  m i(£)  +  irri2(£,)  also  has  the  truncation  property.  The  basic  observation 
is  that  one  can  use  Fourier  series  to  write  any  symbol  of  the  form 

/mi(f)  +wn2(fh 
U - I - ) 
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as 

j,k(^L 

where  ipj{x)  :=  e2mixip  and  ip  is  some  bump  function  which  equals  1  on  the  one-dimensional 
projections  of  the  support  of  ip.  Since  the  Ce  norm  of  ipj  grows  polynomially  in  j  and  ip(j,  k) 
decays  rapidly  in  j,  k  (if  ip  is  sufficiently  smooth),  we  are  done  since  the  product  of  two  LP 
multipliers  is  still  an  Lp  multiplier. 

3.  Nonlinear  hyperbolic  conservation  laws 

Having  developed  our  averaging  lemmata,  we  now  present  some  applications  to  nonlinear 
PDE.  We  begin  with  the  study  (real-valued)  solutions  pit,  x )  =  p(t,  xi, . . .  ,  xf)  G  Z/°°((0,  oo)  x 
Rj()  of  multidimensional  scalar  conservation  laws 

0  d  ^ 

(3.1)  —p(t,x)  +  ^2—Aj(p(t,x))  =  0,  in  Z>'((0,oo)  x  R*). 

3= 1  j 

We  abbreviate  (3.1)  as  pt  +  Vx  •  A(p)  =0  where  A  is  the  vector  of  C'2,e-spatial  fluxes,  A  := 
(Ai,  A2, . . .  ,  Ad). 

Let  x7(v)  denote  the  velocity  indicator  function 

(  1  if  0  <  v  <  7 

X'yi.v)  =  <  -1  if  7  <  V  <  0 
(  0  otherwise. 

We  say  that  p(t,x )  is  a  kinetic  solution  of  the  conservation  law  (3.1)  if  the  corresponding 
distribution  function,  Xp(t,x)  (v),  satisfies  the  transport  equation 

(3.2)  dtXp(t,x){v)  +  a(v)  ■  VxxP(t,x)(v)  =  dvm(t ,  x,  v)  in  V((0,  oo)  x  Rf  x  Mt)), 

for  some  nonnegative  measure,  m(t,x,v )  G  Al  +  ((0,cx))  x  x  M„).  Here,  a(u)  is  the  vector 
of  transport  velocities,  a(u)  :=  (ai(u),...  ,aa(v ))  where  aj(-)  :=  Aj(-),j  =  1,2 ,...d.  The 
regularizing  effect  associated  with  the  proper  notion  of  nonlinearity  of  the  conservation  law  (3.1) 
was  explored  in  [LPT94a]  through  the  averaging  properties  of  an  underlying  kinetic  formulation. 
For  completeness,  we  include  here  a  brief  description  which  will  serve  our  discussion  on  nonlinear 
parabolic  and  elliptic  equations  in  the  next  sections  and  we  refer  to  [LPT94a]  for  a  complete 
discussion. 

The  starting  point  are  the  entropy  inequalities  associated  with  (3.1), 

dtr](p{t,  x))  +  Vx  •  A v(p(t,  x))  <  0  in  V'((0,  oo)  x  R^). 

Here,  rj  is  an  arbitrary  entropy  function  (a  convex  function  from  R  to  R)  and  Av  :=  (Aj?, . . .  ,  Afff) 
is  the  corresponding  vector  of  entropy  fluxes,  Aj(p)  :=  fp  r/'(s)Aj(s)ds,  j  —  1,2, ...  d.  A 
function  p  G  L°°  is  an  entropy  solution  if  it  satishes  the  entropy  inequalities  for  all  pairs 
(77,  An)  induced  by  convex  entropies  rj.  Entropy  solutions  are  precisely  those  solutions  which  are 
realizable  as  vanishing  viscosity  limit  solutions  and  are  uniquely  determined  by  their  L°°  fl  Id- 
initial  data,  po(x),  prescribed  at  t  —  0,  e.g.,  [Lax73].  A  decisive  role  is  played  by  the  one- 
parameter  family  of  Kruzkov  entropy  pairs,  (rj(p;v),A n(p-,v)),  parameterized  by  v  G  R, 

V(p\ v )  :=  I P  ~  v\,  AVj  (p;  v)  :=  sgn (p  -  v)(Aj(p)  -  Aj(v)). 
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Kruzkov  entropy  pairs  lead  to  a  complete  Id-theory  of  existence,  uniqueness  and  stability  of 
first-order  quasilinear  conservation  laws,  [Kr70].  We  turn  to  the  kinetic  formulation.  We  define 
the  distribution,  m(t,x,v )  =  mp(t,x)(v)  by  the  formula 


(3.3) 


m(t,  x,  v)  :=  — 


d, 


v(p’, v )  ~  h(0;  v) 


+  va 


Affp-,v)  -  A'\0;v) 


The  entropy  inequalities  tell  us  that  the  distribution  m  =  mp  is  in  fact  a  nonnegative  measure, 
m{t,x,v )  G  Al+((0,oo)  x  Rf  x  R„).  Next,  we  differentiate  (3.3)  with  respect  to  v:  a  straight¬ 
forward  computation  yields  that  Xp(t,x)(v )  satisfies  the  kinetic  transport  equation  (3.2).  This 
reveals  the  interplay  between  Kruzkov  entropy  inequalities  and  the  underlying  kinetic  formu¬ 
lation;  for  nonlinear  conservation  laws,  kinetic  solutions  coincide  with  the  entropy  solutions, 
[PT91,  LPT94a].  Observe  that  by  velocity  averaging  we  recover  the  macroscopic  quantities 
associated  with  the  entropy  solution  p, 


XP{v)4>{v)dv  =  <h(p) 

where  4>(p)  :=  fp=0  <f>(s)ds  is  the  primitive  of  q b.  In  particular  one  can  recover  p  itself  by  setting 
(p(v)  =  l[_MlM](.y,  M  =  \\p\\L°°- 

We  now  use  the  averaging  lemma  (2.2)  to  study  the  regularity  of  p.  To  this  end  we  first 
extend  (3.2)  over  the  full  Rt  x  R^-space,  using  a  Cq>( 0,  oo)-cut-off  function,  ij}  =  1  for  t  >  e,  so 
that  f(t,x,v )  :=  xP(t,x){v)^{t)  and  g{t,x,v)  :=  mp\t,x,v)il)(t)  +  Xp(t,x)(v)dtff(t)  satisfy 

d tf{t,  x,  v )  +  a(u)  •  Vxf(t,  x,  v )  =  dvg(t,  x,  v )  in  T>'(M.t  x  R)  x  R.t,),  g  G  A4(Rt  xljx  R„). 

Set  /  :=  [inf  p0,  sup  p0]  and  assume  that  the  first-order  symbol  is  non-degenerate,  (2.4),  (2.14), 
namely, 

(3.4)  3a  G  (0, 1)  s.t.  sup  jna(£;  <5)|  <  5ai:  fla(£;  5)  :=  \v  G  /  :  |r  +  a(u)  •  £|  <  <f), 

r2+|5|2=l  1  J 

and 


(3.5) 


3 p  G  [0, 1]  s.t.  sup  sup  |a'(n)  •  £|  <  (W 
|5|=1  Qa(5;<5) 


We  apply  the  averaging  lemma  2.2  for  first-order  symbols,  k  —  1,  with  q  —  l,p  —  2  and 

s’r{%  x  Rf), 


a  =  0,  to  find  that  f(x)  and  hence  p(x)  belong  to  Wloc 
P(t ,  x)  G  Wj£((e,  oo)  x  R)),  s<6a:= 


a 


r  := 


ol  T  4  —  2p 


a  T  4  —  2  p  ol  T  2  —  p 

At  this  stage,  we  invoke  the  monotonicity  property  of  entropy  solutions,  which  implies  that 
for  s  <  1,  ||p(£,  •)llws,1(Rd)  nonincreasing,  and  we  deduce  that  p(t,  ■)  G  kTs,1(R'))  for  t  >  e. 
We  conclude  that  the  entropy  solution  operator  associated  with  the  nonlinear  conservation  law 
(3.1), (3. 4),  pq(-)  i — >  p(t,  •),  has  a  regularizing  effect,  mapping  L°°(Rf)  into  IT)) 


rs,l( 
'  loc 


tg). 


vt  >  e  >  0  :  p(t, .)  e  Wg(Ri),  s  <  »„  Sl  :=  = 


a 


ol  T  4  —  2  p 


Next,  we  use  the  bootstrap  argument  of  [LPT94a,  §3]  to  deduce  an  improved  regularizing 
effect.  The  IT)*’)  (Rt  x  R))-regularity  of  p(t,  x)ijj(t)  implies  that  f(t,  x,  v)  =  Xp(t,x)  (n)V’(t)  belongs 
to  L1(WS,1( Rt  x  Rf),R„);  moreover,  since  dvxP{v)  is  a  bounded  measure,  /  G  L1(hTs,1(Rt  x 
R^),R„)  D  Zd(R t  x  Rf,  hTs,1(Ri,))  and  hence 

/  G  Wj£( Rt  x  Rdx  x  R„),  Vs  <  Sl. 


16 


EITAN  TADMOR  AND  TERENCE  TAO 


Interpolation  with  the  obvious  L°°-bound  of  /  then  yields  that  /  G  Wfj/i Rt  xKjx  M,,)  for  all 
s  <  s i/2.  Therefore,  the  averaging  lemma  2.2  applies  to  /  =  with  q  —  l,p  —  2  and 

a  =  s i/2,  implying  that  pit,  •)  has  improved  W'^! -regularity  of  order  s  <  s2  =  (l  —  9a)si/2  +  9a. 
Reiterating  this  argument  yields  the  fixed  point  ,sy  |  soo  =  20Q/  (1  +  0a)  and  we  conclude  with 
a  regularizing  effect 

Vf  >  e  >  0  :  p0  G  L°°  n  L1^)  ^  pit,  •)  G  5  <  — % - . 

CM  ~r  Z  —  fl 

As  indicated  earlier  in  remark  2.7,  in  the  generic  case,  p  —  1  —  a, 

(3.6)  sup  sup  |a'(u)  •  £|  <  h1_a, 

|5|  =  1  {i)S/:|T+a(?;)-5|<<5} 

which  yields  IVs’ 1  -regularizing  effect  of  order  s  <  a/{ 2a  +  1).  This  improves  the  previous 
regularity  result  [LPT94a,  Theorem  4]  of  order  s  <  a /{a  +  2),  corresponding  to  p,  —  0.  We 
can  extend  the  last  statement  for  general  Lploc  initial  data.  Recall  that  the  entropy  solution 
operator  associated  with  (3.1)  is  Id-contractive.  We  now  invoke  a  general  nonlinear  inter¬ 
polation  argument  of  J.-L.  Lions,  e.g.,  [TaOO,  Interpolation  Theory,  Lecture  8];  namely,  if  a 
possibly  nonlinear  T  is  Lipschitz  on  X  with  a  Lipschitz  constant  Lx  and  maps  boundedly 
Y\  i — >  Y2  with  a  bound  By,  then  one  verifies  that  the  corresponding  K -functionals  satisfy 
K(Tx,t]  X,Y2)  <  LxK(x,  tBy/Lx',  X,  Yi),  and  hence  T  maps  [X,Y-i\e,q  ^  [X,Y2\0,q.  Conse¬ 
quently,  the  entropy  solution  operator  maps  [Ll,  L°°]gy(]  i— >  [Ll,  WfZ \e,q,  0  <  6  <  1  <  q,  and  we 
conclude 

Corollary  3.1.  Consider  the  nonlinear  conservation  law  (3.1)  subject  to  Lp  fl  L1 -initial  data, 
p(0,x)  =  po(x).  Assume  the  non-degeneracy  condition  of  order  a,  (3.4),  (3.6)  holds  over 
arbitrary  finite  intervals  I.  Then  p(t,x )  gains  a  regularity  of  order  s/p' , 

Vt  >  e  >  0  :  po  G  Lp  fl  L\Rdx)  ->  p{t,  •)  G  W£(R^),  a  <  “  1}  ■ 

The  study  of  regularizing  effects  in  one-  and  two-dimensional  nonlinear  conservation  laws  has 
been  studied  by  a  variety  of  different  approaches;  an  incomplete  list  of  references  includes 
[0163] ,  [Ta79] ,  [Ta87] ,  [EE93] ,  [TRB05] . 

We  close  this  section  with  three  examples.  Let  l  >  1  and  consider  the  one-dimensional 
conservation  law 

(3.7)  =  °»  Po  e  [— M, M\. 

It  satisfies  the  non-degeneracy  condition  (3.4)  with  a  =  l/i,  hence  pit,  -)|t>e  G  with 

s  <  a /{2a  +  1)  =  1/(1  +  2).  It  is  well-known,  however,  that  the  entropy  solution  operator  of 
the  inviscid  Burgers’  equation  corresponding  to  t  —  1,  maps  L°°  >  BV,  [0163].  This  shows 

that  the  regularizing  effect  of  order  a/ (2a  +  1)  stated  in  corollary  3.1  is  not  sharp  (although 
the  averaging  argument  is!  consult  [DLW05]  following  [JP02]).  Accordingly,  it  was  conjectured 
in  [LPT94a]  that  (3.4)  yields  a  regularizing  effect  of  order  a. 


Next,  let  i,  m  >  1  and  consider  the  two-dimensional  conservation  law 


d  d 

7 TTP{t,x)  +  — 


dx\  1 1  +  1 


W<+1(i,x)}  +  A 


dx2  lm+1 


pm+1(t,x)  =0,  po  G  [-M,  M\. 
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If  d  m  then  (3.4)  is  satisfied  with  a  =  min{|,  and  we  conclude  pit  >  e,  ■)  E  W^ffL1)  with 
s  <  minl^,  ^b_}.  If  i  =  m,  however,  then  there  is  no  regularizing  effect  since  t' +veff1+vmff2  = 
0  for  r  =  0,^i  +  £2  =  0;  indeed,  pffx  —  y)  are  steady  solutions  which  allow  oscillations  to  persist 
along  x  —  y  =  const.  Other  cases  can  be  worked  out  based  on  their  polynomial  degeneracy;  for 
example, 

d  c)  d  1 

o^p(^x) +  ^rsin(p(M))  +  =  °>  p°e 

has  a  non-degeneracy  of  order  a  =  1/4,  yielding  Ih/// -regularity  of  order  s  <  1/6. 

4.  Nonlinear  degenerate  parabolic  equations 

We  are  concerned  with  second-order,  possibly  degenerate  parabolic  equations  in  conservative 
form 


(4.1) 


d_ 

dt 


d 

■p(t,  x)  + 

3= 1 


() 

dxj 


Aj(p(t,x )) 


a  d2 

a  n  Bjk(p(t,  x))  =  0,  in  V'((0,  00)  x  K*). 

j,k=  1  °  j  k 


We  abbreviate,  pt  +  Vx- A(p)  +  trace  (VX®VX  B(p))  =  0  where  B  is  the  matrix  B  :=  k_1. 

By  degenerate  parabolicity  we  mean  that  the  matrix  B'(-)  is  non-negative,  (B'(-)£,£)  >  0,  E 
Rd.  Our  starting  point  are  the  entropy  inequalities  associated  with  (4.1),  such  that  for  all  convex 
r/s, 

(4.2)  dtrj(p(t,x))  +  Vx  •  A1(p(t,x))  -  trace ^Vx  ®  Vx  Bv(p(t,x))^  <  0  in  D'((0,  00)  x  if ). 

Here,  A'1  is  the  same  vector  of  hyperbolic  entropy  fluxes  we  had  before,  Av  =  (Aff  . . .  ,  A^)  and 
B'  is  the  matrix  of  parabolic  entropy  fluxes,  B''  :=  (Bjk)dk=1,  Bdk(p)  :=  Jp  ff  {s)B'-k{s)ds.  We 
turn  to  the  kinetic  formulation.  Utilizing  the  Kruzkov  entropies,  rj (p',v)  :=  \p  —  v\,  we  define 
the  distribution,  m(t,x,v )  =  (v) , 


m(t,x,v )  := 


g  rj{p-,v)-rj(0]v)  (  An(p;  v)  -  A^O;  v) 

t  X 


(4.3) 


/  „  Bffp;v)  -B^(0;n)\ 

+  trace  (8)  Vx - - - J. 


The  entropy  inequalities  tell  us  that  m(t,x,v )  G  Al  +  ((0,oo)  x  x  R„)  and  differentiation 
with  respect  to  v  yields  the  kinetic  formulation, 


(4.4)  dtxp(t,x){v)  +  a(u)  •  VxxP(t,x)(v)  -  Vj  ■  b (v)VxxP(t,x)H  =  dvm(t,  x,  v ), 

for  some  nonnegative  m  E  A4+  which  measures  entropy -(-dissipation  production.  Here,  a 
is  the  same  vector  of  velocities  we  had  before,  a  =  A',  and  b  is  the  non-negative  diffusion 
matrix,  b  :=  B'  >  0.  The  representation  rj(p)  —  7/(0)  =  f  ff(s)xP(s)ds  shows  that  the  kinetic 
formulation  (4.4)  is  in  fact  the  equivalent  dual  statement  of  the  entropy  inequalities  (4.2).  But 
neither  of  these  statements  settles  the  question  of  uniqueness,  except  for  certain  special  cases, 
such  as  the  isotropic  diffusion,  Bjk(p)  =  B(p)5jk  >  0,  e.g.,  [Ca99],  or  special  cases  with  mild 
singularities,  e.g.,  a  porous-  media  type  one-point  degeneracy,  [DiB93,  Ta97].  The  extension 
of  Kruzkov  theory  to  the  present  context  of  general  parabolic  equations  with  possibly  non¬ 
isotropic  diffusion  was  completed  only  recently  in  [CP03],  after  the  pioneering  work  [VH69]. 
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Observe  that  the  entropy  production  measure,  m,  consists  of  contributions  from  the  hyperbolic 
entropy  dissipation  and  the  parabolic  dissipation  of  the  equation,  m  =  tua  +  iRb-  The  solutions 
sought  by  Chen  and  Perthame  in  [CP03],  p  G  L°°,  require  that  their  corresponding  distribution 
function  \p  satisfies  (4.4)  with  a  restricted  form  of  parabolic  defect  measure  m&\  the  restriction 
imposed  on  me  reflects  a  certain  renormalization  property  of  the  mixed  derivatives  of  p  (or  more 
precisely,  the  primitive  of  \Jh(p)).  Accordingly,  we  can  refer  to  these  Chen-Perthame  solutions 
as  renormalized  solutions  with  a  kinetic  formulation  (4.4).  These  renormalized  kinetic  solutions 
admit  an  equivalent  interpretation  as  entropy  solutions,  [CP03]  and  as  dissipative  solutions, 
[PS05].  A  general  L1-theory  of  existence,  uniqueness  and  stability  can  be  found  in  [CK04],  For 
a  recent  overview  with  a  more  complete  list  of  references  on  such  convection-diffusion  equations 
in  divergence  form  we  refer  to  [Ch03] .  The  regularizing  effect  of  such  equations,  however,  is  less 
understood.  In  [LPT94a,  §5]  we  used  the  kinetic  formulation  (4.4)  to  prove  that  the  solution 
operator,  p0  l— >  p{t,  •)  is  relatively  compact  under  a  generic  non- degeneracy  condition 

sup  |n£(£;0)|  =  0,  fi£(£;0):={u  :  t  +  a(v)  ■  £  =  0,  (b(u)f,f)  =  0}  . 

ki=i 

A  general  compactness  result  in  this  direction  can  be  can  be  found  [Ge90] .  We  turn  to  quantify 
the  regularizing  effect  associated  with  such  kinetic  solutions.  We  emphasize  that  our  regularity 
results  are  based  on  the  ‘generic’  kinetic  formulation  (4.4),  but  otherwise,  they  are  independent 
of  the  additional  information  on  the  renormalized  Chen-Perthame  solutions  encoded  in  their 
entropy  production  measure  m.  The  extra  restrictions  of  the  latter  will  likely  to  yield  even 
better  regularity  results  than  those  stated  below.  We  divide  our  discussion  into  two  stages, 
in  order  to  highlight  different  aspects  of  degenerate  diffusion,  in  Section  4.1,  and  the  coupling 
with  nonlinear  convection,  in  Section  4.2. 

4.1.  Non-isotropic  degenerate  diffusion.  We  consider  the  parabolic  equation 

d  d  d2 

(4.5)  —pit,  x)~Y^  Qx  .dxkBjk^  =  °’  hl  V'^0,  °°)  X 

j,k=  1  -1 

Here,  we  ignore  the  hyperbolic  part  and  focus  on  the  effect  of  non-isotropic  diffusion.  The 
corresponding  kinetic  formulation  (4.4)  extended  to  the  full  Kj  x  x  R„  reads 

d tf(t,  x,  v )  -  Vj  •  b (v)Vxf(t,  x,  v )  =  dvg(t,  x,  v ),  /  :=  yy0(f),  g  G  M+(Rt  x  R!^  x  R„). 

Set  /  :=  [inf  p0,  sup  po]-  The  corresponding  symbol  is  C{r,i£,v)  =  ir  +  (b(u)£,£)  and  it  suffices 
to  make  the  non-degeneracy  assumption  (2.4),  on  the  second-order  homogeneous  part  of  the 
symbol  £(0,  i£,v)  =  (b(u)£,  £).  We  make 


(4.6)  3a  G  (0, 1)  s.t.  sup  |flb(£;  <5)1  ^  Sa,  (£;  5)  :=  {u  G  /  :  0  <  (b(u)^,^)  <  5}. 

id=i 

and 

(4.7)  3p  G  [0, 1]  s.t.  sup  sup  |(b,(u)£,  £)|  < 

|5|=1  Ob(5;5) 


We  apply  the  averaging  result  2.3  with  q  =  l,p  =  2,  a  =  0  and  k  = 
hence  p(t,x)  belong  to  W^(Rt  x  Rf), 

P(t,  •)  G  W£((e,  oo)  x  Rf),  8  <  2 9a,  6a  :=  — -, 


2,  to  find  that  f(t,x )  and 


o  1  —  2p 

(\  2  —  p, 


r  := 


VELOCITY  AVERAGING  AND  REGULARITY  OF  QUASILINEAR  PDES 


19 


We  follow  the  hyperbolic  arguments.  The  kinetic  solution  operator  associated  with  (4.1)  is 
Id-contractive,  hence  || pit,  •)llws,1(Rd)  nonincreasing  and  we  conclude  that  Vt  >  e,  p(t,  •)  has 

-regularity  of  order  s  <  si,  Si  :=  2 9a  =  2a /(a  +  4  —  2/i).  We  then  bootstrap.  Since 
Xp(v)^(t)  G  W/j/  (Rt  x  x  R„)  for  all  s  <  s i/2,  we  can  apply  the  averaging  lemma  2.1 
with  a  =  si/2  leading  to  IT^1 -regularity  of  order  s2  (1  —  9a)s\/2  +  2 9a  with  fixed  point 

sk  T  ^oo  =  4 9a/ (1  +  9a), 

(4.8)  Vi>e>  0:  »  <  — d — ■ 

a  +  2  —  p, 

We  distinguish  between  two  different  types  of  degenerate  parabolicity,  summarized  in  the 
following  two  corollaries. 

Corollary  4.1.  [Degenerate  parabolicity  I.  The  case  of  a  full  rank]  Consider  the  degenerate 
parabolic  equation  (f-1)  subject  to  L°°  CLl-initial  data,  p(0,x)  =  p0-  Let  Ai(v)  >  A2(u)  > 
...A  d(v)  >  0  be  the  eigenvalues  of  b(v)  and  assume  that  A  d(v)  ^0  over  I  =  [inf  p0,  supp0]- 
Then,  >  Arf(u)|^|2  and  p(t,x )  has  a  regularizing  of  order  s  <  2a/ {a  +  2  —  p)p' ,  i.e., 

(4-8)  holds  with  a  and  p,  dictated  by  \<i(v)  =  Arf(b(u))? 

|Oa(<5) |  <  and  sup  sup  ((b^u)^ ^)|  <  (P,  Q\(5)  \=  {v  E  I  :  0  <  \d(v)  <  5}. 

|^|=i 

Corollary  4.1  applies  to  the  special  case  of  isotropic  diffusion, 

(4.9)  ^p(t,  x)  -  A B(p(t,  x))  =  0,  B\v)  >  0, 

subject  to  L°°  fi  L1 -initial  data,  p(0,  •)  =  po-  If  b(-)  :=  B'{v)  is  degenerate  of  order  a  in  the 
sense  that, 

|Ofo(<5) |  <  Sa,  and  sup  \br (v)\  <  d1-01 ,  Llb(8)  :=  {v  G  I  :  0  <  b(v)  <  5}, 
venb(5) 

then  corollary  4.1  implies  Vt  >  e  :  pit ,  •)  G  Wfj],  s  <  2a/i/2a  +  1).  For  Lp  fl  L1-data  p0,  the 
corresponding  solution  pit,  ■)  gains  My (f' -regularity  of  order  s  <  2a /(2a  +  l)p'  and  we  conjec¬ 
ture,  in  analogy  with  the  hyperbolic  case,  that  the  non-degeneracy  (4.6)  yields  an  improved 
regularizing  effect  of  order  2 a/p' .  Existence,  uniqueness  and  regularizing  effects  of  the  isotropic 
equation  (4.9)  were  studied  earlier  in  [BC79,  BC81a,  BC81b].  The  prototype  is  provided  by 
the  porous  media  equation, 

d  1 

(4-10)  — p(t,x)  -  A{-^-y|pn(t,x)|p(t,x)}  =  0,  p(0,  x)  =  po(x)  >0,  po  G  L°° . 

The  velocity  averaging  yields  M-7^1 -regularity  of  order  2/(n  +  2)  and,  as  in  the  hyperbolic  case, 
it  does  not  recover  the  optimal  Holder  continuity  in  this  case,  e.g.  [DiB93,  Ta96].  In  fact, 
the  kinetic  arguments  do  not  yield  continuity.  Instead,  our  main  contribution  here  is  to  the 
non-isotropic  case  where  we  conjecture  the  same  gain  of  regularity  driven  by  A^(b(u)),  as  the 
isotropic  regularity  driven  by  b(v). 

We  continue  with  the  more  subtle  case  where  b(-)  does  not  have  a  full-rank,  so  that 

31,  1  <  i  <  d  :  Ai(u)  >  . . .  \e(v)  >  0,  Am(u)  =  . . .  =  Xd(v)  =  0 

Despite  this  stronger  degeneracy,  there  is  still  some  regularity  that  can  be  ‘saved’.  To  demon¬ 
strate  our  point,  we  consider  the  2D  case. 
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Corollary  4.2.  [Degenerate  parabolicity  II.  The  case  of  a  partial  rank]  We  consider  the  2D 
degenerate  equation 

d  r  d2  d2  d2  'i 

(4-n)  mp{t’x)  ~  \d^BlMtiX))  +  d^2BMx))  +  dAB^x)))  =  °- 

subject  to  L°°  fi  L1  -initial  data,  p(0,  x)  =  p0.  Assume  strong  degeneracy,  b\2(v )  =  4:bu(v)b22(v) , 

so  that  \2(h(v))  =  0,  \/v  G  /  =  [inf  p0,  sup0].  In  this  case,  (b(i>)£,£)  =  ^v/5n(v)^i  +  v/^22('t;)^2^ 

and  p(t,  •)  admits  a  Wfj] -regularity  of  of  order  s  <  2a/ {ex  +  2  —  p)  which  is  dictated  by  the  non¬ 
degeneracy, 

(4.12)  sup  nb(f;5)  ZSa,  :=  \v  E  I  :  |\/&ii(u)fi  +  y/b^iv)^  <  <$}, 

l€l=i  L  J 

and 


(4.13) 


sup  sup 
l€l=i  venh(t-,6) 


Vll(v)  c  b22^)  t 

Vbu(v)  1  Vh2(v)  2 


<  5»~1/2. 


We  distinguish  between  two  extreme  scenarios. 

(i)  If  |6n(u)|  S>  | b22 (v) |  \/v  G  I,  then  the  regularizing  effect  (4.8)  holds  with  (a,  p)  dictated 
by  b22(v), 

|^&22(<S)|  ~  ba  and  sup  \b22(v)\  <  5^,  klb22(S)  (u  G  /  :  0  <  b22(v)  <  £}• 

ueO{,22  (6) 


(ii)  If  bu(v)  =  b22(v)  \/v  G  I,  then  there  is  no  regularizing  effect  since  the  symbol  ^v/5n(u)£i  + 
\  2 

Vbf/iv)^)  vanishes  for  all  £1  ±  £2  =  0  (so  that  (4.12)  is  fulfilled  with  a  =  0).  Indeed,  the 
equation  (4.11),  with  &n(u)  =  b22{v)  =:  B\v),  takes  the  form 


0, 


and  we  observe  that  p0(x  =F  y)  are  steady  solutions  which  allow  for  oscillations  to  persist  along 
x  =F  V  =  const. 


4.2.  Convection-diffusion  equations.  We  begin  with  the  one-dimensional  case 
(4.14)  ^p{t,  x )  +  ^-A(p(t,  x))  -  ^B(p(t,  x))  =  0. 

We  consider  the  prototype  example  of  high-order  Burgers’  type  nonlinearity,  a(v)  :=  ve,£  >  1 
combined  with  porous  medium  diffusion  b(v)  =  |u|",n  >  1.  The  corresponding  symbol  is  given 
by  £((r,  *£),  v)  =  ir+vHf  +  |u|"£2.  We  study  the  regularity  of  this  convection-diffusion  equation 
using  the  averaging  lemma  2.3,  which  employs  the  size  of  the  set 

nc(j ;  S)  :=  {v  I  J\t  +  vef\  +  J2\v\n£2  <  <j}  ,  r2  +  f2  =  1,  J  >  1,  5  <  1. 

Comparing  diffusion  vs.  nonlinear  convection  effects,  we  can  distinguish  here  between  three 
different  cases.  Clearly,  Dc(J]S)  C  Db  :=  {v  :  \v\n  <  S/J2},  hence  uc(J;6)  <  {5/ J2)l^n  and 
(2.19)  holds  with  ab  =  1  /n  and  /3b  =  2.  We  shall  use  this  bound  whenever  n  <  I,  which  is  the 
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case  dominated  by  the  parabolic  part  of  (4.14).  Indeed,  in  this  case  we  have  {8  /  J)1^  >  {5/  J2)l^n 
which  in  turn  yields 

sup  |£v((t,*£),u)|  £  sup  (J\v\e~l  +  J2\v\n~l)  <  J2/n51_1/n. 

v&Cl  v&Vlb 

This  shows  that  (2.20)  holds  with  (A b,  Pb)  —  (ab,  1  —  otb)  and  velocity  averaging  implies  Wff/- 
regularizing  effect  with  Sobolev  exponent  s  <  /3bab/{3ab  +  2), 

2 


We  also  have,  f C  fla  ■—  {v  :  \v\e  <  8/J},  so  that  uffJ',8)  <  {8/J)1^,  i.e.,  (2.19)  holds 
with  aa  =  l/£  and  (3a  =  1.  We  shall  use  this  bound  whenever  n  >  2d,  which  is  the  case  driven 
by  the  hyperbolic  part  (4.14).  In  this  case,  {8/ J)1^  <  {8/ J2)1/”,  hence 

sup  |£„((r, i£),u)|  <  sup  (J|u|^_1  +  J2|u|n_1)  < 

v£flc(J;8)  v£  fla 


implying  that  (2.20)  is  fulhlled  with  (A a,  (J,a)  =  {aa,  1  —  aa).  The  corresponding  Sobolev  expo¬ 
nent,  s  <  /3aaa/{3aa  +  2),  is  then  given  by, 

1 


Finally,  for  intermediate  n’s,  £  <  n  <  2d,  we  interpolate  the  previous  two  cu^-bounds  (which 
are  valid  for  all  n’s), 

uc{J;  8)<{8/J){-1-W\8/J2)Un, 


for  some  Q  E  [0, 1],  which  we  choose  as  £  :=  (n/£)  —  1,  so  that  (2.19)  holds  with  a  —  (1  —  (/)/£  + 
(t/n,(3a  =  (1  —  ()/£  +  2(t/n  and  (2.20)  holds  with  (A ,  p)  =  (a,  1  —  a).  This  then  yields  the 
Sobolev-regularity  exponent,  s  <  (3a /{3a  +  2), 


n  +  {2d  —  n)( 

3 n  +  3 {£  -  n)C  +  2nd’ 


c  :=  (n/£)  -  1,  £  <  n  <  2d. 


An  additional  bootstrap  argument  improves  this  Sobolev  exponent,  s  <  /3a /{2a  +  1),  and  we 
summarize  the  three  different  cases  in 


Corollary  4.3.  The  convection- diffusion  equation 


r)  f)  r  1  'l  cP  (  1  'l 

(4-15)  fop&x)  +  d^{lTipe+1^x^}  ~  =  °>  po  G 


has  a  regularizing  effect,  p0  E  L°°{M.X)  i— >  p{t  >  e,  ■)  E  Wl 

n  +  {2d  -  n)0,n 


£(R») 


S£,n  — 


2  n  +  2{£  —  n)Qjn  +  n£' 


,  of  order  s  <  spn  given  by 

0  n  <  £ 

Ct,n  ■■=  {  {n/£)  -  1  £<n<2£ 

n  >  2d 


We  note  that  when  n  <  £,  then  Qb  C  £la  and  (4.15)  is  dominated  by  degenerate  diffusion  with  a 
regularizing  effect  of  order  syn  =  ,sn  =  2/{n  +  2).  Thus,  we  recover  the  same  order  of  regularity 
we  met  with  the  ‘purely  diffusive’  porous  medium  equation  (4.10).  If  n  >  2d,  however,  then 
Qa  C  and  it  is  the  hyperbolic  part  which  dominates  diffusion,  driving  the  overall  regularizing 
effect  of  (4.15)  with  order  sy.n  =  se  =  !/{£+ 2);  we  recover  regularity  with  the  same  order  which 
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we  met  with  the  ‘purely  convective’  hyper-Burgers’  equation  (3.7).  Finally,  in  the  intermediate 
‘mixed  cases’,  £  <  n  <  27,  we  find  a  regularity  of  a  (non-optimal)  order 


S£,n 


n£  +  (27  —  n)(n  —  £) 

2 n£  —  2{n  —  £)2  +  nl2  ’ 


We  turn  to  the  multi-dimensional  case  (4.1).  The  regularizing  effect  is  determined  by  the 
size  of  the  set 

n(J;6)  =  {v  |  J|r  +  a(u)-e|  +  J2(b(u)e,0<5},  r2  +  |^|2  =  1,  J  >  1,  6  <  1. 

Assume  that  the  degenerate  parabolic  part  of  the  equation  has  a  full-rank,  so  that  the  smallest 
eigenvalue  of  b(u),  X(v)  =  A<*(b(r>))  satishes 

(4.16) 

|^b(A) |  <  A"b  and  sup  sup  |(b'(u)£,£)|  <  <51_ab,  flb(5)  :=  |u  0  <  Arf(b(u))  <  8 
|g|=ivenb(<5)  1 


In  this  case,  u>c(J',  5)  ~  (5/  J2)ab  which  yields  a  gain  of  Wjf^- -regularity  of  order  s  <  2ab/(2o:b  + 
1).  If  in  addition,  the  hyperbolic  part  of  the  equation  has  a  non-degeneracy  of  order  aa,  namely 

sup  |f2a((r,  £),  5)  |  <  5Qa  and  sup  sup  |a'(u)  •  £|  < 
r2+|5|2=l  151=1  «e!la  (5) 


where  fia((T, £);  h)  :=  |u  :  |r  +  a(u)  •  £|  <  hj,  then  we  can  argue  along  the  lines  of  corollary 

4.3  to  conclude  that  there  is  an  overall  IT))^1 -regularity  of  order  dictated  by  the  relative  size  of 
2ab/(2a:b  +  1)  and  aa/(2aa  +  1).  As  an  example,  we  have  the  following. 


Corollary  4.4.  Consider  the  two-dimensional  convection- diffusion  equation 
(4.17) 


d_ 

dt' 


p(t'x)  +  -k{ihpW[t'x)} +  s^{^rTTp”+1(i,l)}  -  £  d^rf>iMt’x))  =  °- 

j,k=  1  J 


d 


d 2 


with  non- degenerate  diffusion,  B'{y )  >  |u|".  Then,  its  renormalized  kinetic  solution  admits  a 
Wff] -regularizing  effect,  p0  G  L°°  i— >•  p(t,  •)  G  hF^’^M2),  of  order 

1  1 


s  <  < 


Slrn  ■=  mill 


Sn.  • 


1  +  2  ’  m  +  2  / 


n  +  2’ 

[  st  ,m,n  e  [st  ,774  5  ^7 1\  5 


if  n  >  2  max (7,  m)  and  £  f  m, 

if  n  <  min(7,  m)  or  £  =  m, 
if  min (£,m)  <  n  <  2max(7,  m). 


Finally,  we  close  this  section  with  a  third  example  of  a  fully-degenerate  equation 


(4.18) 


4p((.x)  +  (A  +  ^_)^p((>:c)) 


d 2  d 2 

dx\  dx\dx2 


+ 


dxl 


B{p(t,  x))  =  0. 


In  this  case,  there  is  a  stronger,  rank-one  parabolic  degeneracy  with  no  regularizing  effect  from 
the  purely  diffusion  part,  since  (b(i>)£,£)  =  0,  Vfi  —  £2  =  0,  and  no  regularizing  effect  from 
the  purely  convection  part  where  a(u)  •  £  =  0,  V£i  +  £2  =  0.  Nevertheless,  the  combined 
convection- diffusion  does  have  a  regularizing  effect  as  demonstrated  in 
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Corollary  4.5.  Consider  the  two-dimensional  convection- diffusion  equation 


d_ 

dtr 


p(t’x)  +  (ly+CKrr rW(t’x)} 


(4.19) 


( 


d 2 


-2- 


d 2 


+ 


d2 


-\pn(t,x)\p(t,x)  \  =  0. 


\dx2  “ dxidxz  '  dx^Jin  +  l' 

For  n  >  2d  it  admits  a  regularizing  effect,  p0  E  L°°  n  L1  ^  p(t ,  -)  E  IF//!  (M/);  of  order 
s  <  6/(2  +  2n  —  £). 


As  before,  the  Sobolev  exponent  computed  here  is  not  necessarily  sharp,  and  as  a  result  no 
gain  of  regularity  is  stated  for  n  <  2d. 

For  proof,  we  use  the  averaging  lemma  2.3  with  the  usual  (p,q)  =  (2,1),  which  yields  a 
Sobolev-regularity  exponent  of  order  s  =  /3a( 2  —  p  —  A)/(ck  +  4  —  2/i),  where  a,/?,  A  and  p 
characterize  the  degeneracy  of  the  symbol  associated  with  (4.19), 


£(CM0»u)  =  Jir  +  +6)  +  J'2\v\n\fi  -6I2,  r2  +  \f\2  =  1. 


We  consider  first  those  £’s,  |£|  =  1  such  that  |£i  — £2|  <  1/5.  Here  we  have  (say)  |£i+£2|  >  1/10, 
so  that  f lc(J]  il)  C  Oa  =  {v  :  \v\l  <  8/  J},  and  (2.19)  holds  with  (aa,  j3a)  =  {l/d,  1);  moreover, 
the  growth  of  Cv  sought  in  (2.20)  is  bounded  by, 


sup  \Cv({r,if),v)\  <  J 

v£flc(J;S) 


5  \  (<-!)/< 

j) 


+  J2 


/  8  \  (n-i)/* 
VJ7 


< 


JV, 


where 

f  A  :=  l/d  and  p  1  —  l/d,  if  n  >  2d, 

|  X  2  —  (n  —  l)/d  and  p  :=  (n  —  1  )/d,  if  n  <  2d. 


In  particular,  if  n  <  2d  then  the  Sobolev  exponent  vanishes  since  2  —  p  —  A  =  0,  and  we 
cannot  deduce  any  regularizing  effect  in  this  case.  If  n  >  2d,  however,  we  compute  the  Sobolev 
exponent  as  before,  s  =  se  =  1/(2 d  +  3).  Next,  we  consider  the  case  when  |£i  —  £2|  >  1/5, 
so  that  Qc{J]8)  C  Hb  =  {v  ■  |u|n  ~  8/J2}  and  (2.19)  holds  with  (a6,/?6)  =  (l/n,2).  As  for 
(2.20),  we  have 


sup  \Cv((r,i^),v)\  <  J 
vGttciJX) 


8  \h-T)/n 

l2) 


+  J2 


8  \  (ji—l)/n 

J2) 


< 


J2X8f 


where 

f  A  =  1/n  and  p  —  1  —  1/n,  if  n  <  d, 

\  A  =  1/2  —  [d  —  l)/n  and  p  =  (d  —  1  )/n,  if  n  >  d. 

In  particular,  if  n  >  2d  we  compute  in  this  case  a  smaller  Sobolev  exponent  s  —  3/ (3  — 2£+4n)  < 
sc-  The  regularity  result  follows  from  the  bootstrap  argument  we  discussed  earlier  which  yields 
the  final  Sobolev  exponent  s  =  2f3a(2  —  p  —  A)/(a  +  2  —  p).  □ 


5.  Nonlinear  degenerate  elliptic  equations 
We  consider  the  nonlinear,  possibly  degenerate  elliptic  equation 
d  d2 

-  E  7T~: r^W1))  =  S in  ®'<r),  b(.)  :=  B'(.)  >  0. 

j  k=l  OX]OXk 


(5.1) 
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augmented  with  proper  boundary  conditions  along  the  C'1,1-boundary  dY.  We  assume  that  the 
nonlinear  source  term,  S(p),  is  further  restricted  so  that  blow-up  is  avoided. 

We  begin  with  formal  manipulations,  multiplying  (5.1)  against  p'(p),  and  ‘differentiating  by 
parts’  to  find 


02 

-  £  ~  '/(P)S(P)  = 

n  h=  1  J  K 


j  =  l  J  j,k=l 


dp  dp 
dxj  dxk 


v'(p)s(p)  = 


=  — trace (8)  V.r  B'?(p)j  +  r]"(p)(b(p)Vxp,\7xp)  -  r]'(p)S(p). 


We  arrive  at  the  entropy  inequalities  associated  with  (5.1),  stating  that  sufficiently  smooth 
solutions  of  (5.1)  satisfy,  for  all  convex  p’s, 

(5.2)  -trace  (v*  0  Vx  B n(p)^j  <  p/(p)S'(p),  Br?(p)  =  | B^k(p)  :=  J  r}'(s)bjk(s)ds} 

So  far,  we  have  not  specified  the  notion  of  solutions  for  (5.1)  since  it  seems  that  relatively  little 
is  known  about  a  general  stability  theory  for  degenerate  equations  such  as  (5.1).  The  difficulty 
lies  with  the  type  of  degeneracy  which  does  not  lend  itself  to  standard  elliptic  regularity  theory, 
because  the  Bjk  s  degenerate  dependence  on  p,  nor  does  it  admit  the  regularity  theory  for 
viscosity  solutions,  e.g.,  [CIL92],  [CC95],  because  of  their  degenerate  dependence  on  p  rather 
than  Vxp.  We  refer  to  the  works  of  Guan  [Gu97],  [Gu02]  who  shows  that  in  certain  cases,  one 
is  able  to  ’’lift”  a  C'1,1-regularity  of  p  into  a  statement  of  C'0O-regularity.  Using  the  existence 
of  smooth  viscosity  solutions  in  the  uniformly  elliptic  case  where  b(v)  >  A  >  0,  (5.2)  could  be 
then  justified  by  the  “vanishing  viscosity  limit”,  forming  a  family  of  regularized  solutions,  pA 
associated  with  bA(s)  :=  b(s)  +  \Icixd  and  letting  A  |  0+.  Next  comes  the  kinetic  formulation 
of  (5.1)  which  takes  the  form 


(5.3)  -Vj  ■  b (v)VxxP(v)  +  S(v)-^xP(v)  =  ^m(x,v)  in  V'{T  x  R„), 

for  some  nonnegative  m  E  Ai+  which  measures  ’’entropy  production”.  Indeed,  for  an  arbitrary 
convex  “entropy”,  p,  the  moments  of  (5.3)  yield 

0  >  —  J  rj"(v)m(x ,  v)dv  =  j  p  (v)-^-m(x,  v)dv  = 

d  r  Q2  /•  d 

=  -  /  v'(v)h3k(v)-K- Q-  Xp(v)dv+  ri\v)S(v)—xP{v)dv  = 

j,k=  1  ’  i  k  J 

d  d2  f  /  \ 

=  -  gx  g —  Bjk(p)  ~  J  (v'(v)S(v))'xP{v)dv  =  —trace  ^ V x  0  Vx  B'?(p)J  -  rj'(p)S(p). 

j,k=l  ^ 

Thus,  the  kinetic  formulation  (5.3)  is  the  dual  statement  for  the  entropy  inequalities  (5.2).  We 
postulate  that  p  is  a  kinetic  solution  of  (5.1)  if  the  corresponding  distribution  function  Xp{x)(v) 
satisfies  (5.3),  and  we  address  the  regularizing  effect  of  such  kinetic  solutions. 
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To  use  the  averaging  lemma,  we  first  extend  (5.3)  over  the  full  x  iC-space.  Let  if}  be 
C“(i?+)-cut-off  function,  i/j(s)  =  1  for  s  >  e,  and  let  ((x)  denote  the  smoothed  distance 
function  to  the  boundary,  ((x)  =  ip(dist(x,  dT)),  then  f(x,v)  :=  Xp(x)(v)C(x)  satisfies,  in 
V'(Rdx  xl„), 


d2 


j,k= 1 


dxjdxk' 

d  ^  d  d 

C( x)m(x ,  v )  +  E  w—  {bjk(v)(Xk(x)xP(v))  +  jr-(bjk(v)(Xj(x)xp(v))  + 

j  k= i  OX3  OXk 


dv 


(5.4) 


S(v)((x)g-Xf(v)  -  Y  bjt(v)CVk(x)x,(v) 

j,k=  1 

d 


=■  Q^gi(x,v)  +  \vxg2(x,v)  +  g3(x,v)  +  g±(x,v). 


Assume  that  (5.1)  is  non-degenerate  in  the  sense  that  there  exists  an  a  G  (0, 1)  such  that 


(5.5) 


Slbffl 


<  5a  and  sup  sup  |(b/(u)£,£)|  <  51  ", 
|£|=it,enb(<5) 


^b(5)  :=  |u  G  /  :  (b(u)f,f)  < 


We  examine  the  contribution  of  each  of  the  four  terms  on  the  right  of  (5.4)  to  the  overall 
tTs>] -regularity  °f  ./)  appealing  to  the  different  averaging  lemmata  term-by-term.  The  first 
term  on  the  right  involves  the  bounded  measure  g±  =  ( m ;  averaging  lemma  2.3  with  the  usual 
(p,  q )  =  (2, 1)  then  yields  that  the  corresponding  average  f  1  has  a  1  -regularity  of  order 
Si  <  26*! ,  9\  =  a/ (3a  +  2).  The  second  term  on  the  RHS  of  (5.4)  involves  the  gradient 
of  the  uniformly  bounded  term  g2  =  ^2bjk(v)(Xk(x)xP(v)]  here  we  can  use  averaging  lemma 
2.1  with  rj  —  q  —  1  to  conclude  that  the  corresponding  average  f2  has  1T)*2C’  -regularity  of 
order  s2  <  92  =  a/ ( a  +  2)  (in  fact,  with  q  —  2  one  concludes  a  better  WS2,2  regularity  of  order 
s 2  <  a/2).  The  remaining  two  terms  on  the  RHS  of  (5.4)  yield  smoother  averages  and  therefore 
they  do  not  affect  the  overall  regularity  dictated  by  the  first  two.  Indeed,  dvxP(v)  and  hence 
g3(x,v)  =  S(v)((x)dvxP(v)  is  a  bounded  measure  and  averaging  lemma  2.1  with  rj  =  N  =  0 
implies  that  the  corresponding  average  f3  belongs  to  the  smaller  Sobolev  space,  IT/)/’ 1  of  order 
s3  <  2 93,  03  —  a/(a  +  2).  Finally,  the  last  term  on  the  right  of  (5.4)  consists  of  the  bounded 
sum,  bjk(v)CXjXk(x)xp',  with  rj  —  N  —  0  and  q  —  2,  the  corresponding  average  f4  has 

a  W£?  regularity  of  order  s4  <  2d4,  d4  =  a/2. 

Next,  we  iterate  the  bootstrap  argument  we  mentioned  earlier  in  the  context  of  hyperbolic 
conservation  laws.  The  first  IT’S1: '-bound  together  with  the  L°°-bound  of  /  imply  a  IT)/',’2-bound 
with  cr 4  =  si/2,  which  in  turn  yields  the  improved  regularity  of  f1  G  IT/jc ,  s  <  (1  —  Q\)(J\  +  2 9\. 
Thus,  for  the  first  term  we  can  iterate  the  improved  regularity,  s4  i— >•  (1  —  0\)s\/ 2  +  29\, 
converging  to  the  same  fixed  point  we  had  in  the  parabolic  case  before,  s4  <  2a/ (2a  +  1). 
The  second  term  requires  a  more  careful  treatment:  as  we  iterate  the  improved  regularity  of 
XP  ^  IT/// ,  we  can  express  the  term  on  the  right  of  (5.4)  as  A Vsg2  with  rjs  1  —  s  and  with  g2 
standing  for  the  sum  of  L 1  -bounded  terms,  g2  =  As  bjk(v)(XkxP(v).  Consequently,  averaging 
lemma  2.1  yields  the  hxed  point  iterations  s2  (1  —  92) s2/ 2  +  (2  —  rqS2)92  with  limiting  regularity 
of  order  s2  <  a.  The  remaining  two  terms  are  smoother  and  do  not  affect  the  overall  regularity: 
a  similar  argument  for  the  third  term  yields  the  hxed  point  iterations,  S3  1— >  (1  —  93)s3/2  +  293 
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with  a  fixed  point  S3  <  2a/(a  + 1),  while  the  fourth  term  remains  in  the  smaller  Sobolev  space 
Wa’2.  We  summarize  with  the  following  statement. 


Corollary  5.1.  Let  p  G  L°°  be  a  kinetic  solution  of  the  nonlinear  elliptic  equation  (5.1)  and 
assume  the  non- degeneracy  condition  (5.5)  holds.  Then  we  have  the  interior  regularity  estimate 
for  all  D  c  T, 


p(x)eWpc(D),  s< 


a,  if  a  <  1/2, 


2  a 

2a  +  1’ 


if  1/2  <  a  <  1. 
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